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Chiral perturbation theory is nowadays a well-established approach to incorporate the chiral con- 
straints from QCD. Nevertheless, for systems involving one baryon, the power counting which 
dictates the chiral order of observables is not as simple and consensual as in the purely mesonic 
case. The heavy baryon approach, which relies on a non-relativistic expansion around the limit of 
infinitely heavy baryon, recovers the usual power counting but destroys some analytic properties 
of the scattering amplitude. Some years ago, Becher and Leutwyler proposed a Lorentz-invariant 
formulation of chiral perturbation theory that maintains the required analytic properties, but at the 
expense of a less intuitive power counting. 

Aware of the shortcomings of the heavy baryon formalism, the Sao Paulo group derived the two- 
pion exchange component of the nucleon-nucleon potential in line with the works of Becher and 
Leutwyler. A striking result was that the long distance properties of the potential is determined by 
the specific low energy region of the pion-nucleon scattering amplitude where the heavy baryon 
expansion fails. In this talk I will discuss the origin of such failure and how it reflects in the 
asymptotics of the nucleon-nucleon interaction. Some results for phase shifts and deuteron prop- 
erties will be shown, followed by a comparison with the heavy baryon predictions. 
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1. Analyticity constraints 

Becher and Leutwyler [jl|] demonstrated that, in the baryon sector, convergence of the chiral 
expansion of certain quantities (form factors, for instance) is a delicate issue. This is closely related 
to the presence of an anomalous threshold in the triangle diagram (l.h.s. of Fig.[l]) for momentum 
transfer t close to the two-pion threshold, 4m|. Such analytic structure is completely missed in the 
usual heavy baryon (HB) formulation of chiral perturbation theory (#PT). It can only be recovered 
via a resummation of the HB series to all orders. To understand this issue (for more details, see for 
instance Ref. [0]) one starts with the spectral representation of the triangle integral, 
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The argument x = Im^ \ft' — 4m 2 /(t' — 2m 2 ) is formally counted as order q~ 1 in the HB expansion, 

which yields tan _1 x = n/2—\/x+\ /3x 3 H .In fact, the first two terms reproduce the HB result 

for the triangle graph, 
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where q = \q\, and L(q) and A(q) are the usual HB loop functions, 



L{q)= VM±l ln VM±l±l, A(,) = larctan-<L. (1.3) 

q 2m K 2q 2m n 

However, it does not take into consideration the case \x\ < 1, where t' gets closer to 4mr n . This re- 
gion controls the long distance behavior of the triangle diagram, as can be seen by its representation 
in configuration space p, H, B], 



r{r) = i r dt >[ ' r di c± Imr(0 . (W) 

nhml J (2ny t' + q 2 \K 2 hm\ 



Clearly one sees that, in order to have a good asymptotic description of T(r), one needs a decent 
representation for Im y(f') near t' = 4m|, which cannot be provided by the usual HB formulation. 

We want to stress that the expansion in l/m^ of our two-pion exchange nucleon-nucleon 
potential (TPEP) should, in principle, recover the expressions from HB^PT. This can be used as a 
cross-check of our calculation (see next section). However one must keep in mind that, due to the 
problem described above, such an expansion is ill-defined at long distances and should be avoided. 

2. The Lorentz-invariant two-pion exchange amplitude 

The two-pion exchange (TPE) amplitude that contributes to the NN interaction is closely 
linked to the 7rN amplitude, which is generically expressed in terms of two invariant amplitudes 
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Figure 1: The triangle diagram that contributes to the nucleon scalar form factor. On the right panel are 
shown ratios of the two first terms of y(— <? 2 )|hb> divided by y(— <? 2 )|bl- F° r more details, see Ref. 



D ± (v,0 and B ± (v,t), 



^iN — Sab + iSbac X c T ffN , 



1 jzN 



u(p') 



D 



2niN 



u(p'), 



(2.1) 



where p and p' are the initial and final momenta of the nucleon, respectively, while v = [(p + 
k) 2 — (p — k') 2 ]/4mN and t = (k' — k) are the usual Mandelstam variables. This allows the TPE 
amplitude to be written as 
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where the superscript (i) refers to nucleon i. The symbol /[•••] represents the four-dimensional 
integration with two pion propagators, 



d 4 Q 1 



1 



(2nY (P-ml) (kP-n%)' 



(2.4) 



with Q= (k + k') /2 the average momentum of the exchanged pions. The profile functions J^'s are 
covariant loop integrals written in terms of the amplitudes D ± and B ± , 



^bd = -if[---]Q x BU±DW±, = -!/[•■ -}Q l Q p BW±BW±. 



(2.5) 



With the 7iN amplitude to 0(q i ) as input [[I]], one generates the TPE amplitude to 0(q 4 ), repre- 
sented graphically by Fig.0. The contributions are grouped in three families of diagrams, according 
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to their topology. The first line of Fig|2| corresponds to the irreducible one loop graphs with vertices 
from the 0(q i ) 7iN chiral Lagrangian, Jtf£$, with coupling constants at their physical values (fam- 
ily I). The second line (family II) contains two-loop diagrams with an intermediate Tin scattering, 
while the third line (family III) comprises one loop graphs with vertices from Jz^ N and Jz^ N , as 
well as one-loop vertex corrections, parametrized in terms of ttN subthreshold coefficients. 

The diagram with a plannar box topology (second of family I) contains a reducible piece 
(iterated one pion), which has to be subtracted in the definition of the potential. This subtraction 
is represented by the third graph of family I, the symbol "+" standing for a nucleon with only its 
positive-energy projection. It is well-known that such projection is not uniquely defined [g]. In 
order to compare our results with the HB ones we adopt the same subtraction followed by Ref. [0]. 




Figure 2: Dynamics of the Lorentz-invariant TPER 



3. Comparison with the HB approach 

The expressions for the Lorentz-invariant TPEP at N 3 LO is lengthy and will not be reproduced 
here (see ||] instead). They are written in terms of Lorentz-invariant loop integrals which, due 
to the reasons mentioned in SecJIJ, do not admit the naive HB expansion. Nevertheless, if one for- 
mally performs such unjustified expansion, one recovers all but three terms of the HB expressions 
[§]. These discrepancies come from two-loop diagrams of family II and affect mostly the central 
isovector component. Numerically, the difference is about 10% up to 3fm and decreases to 5-2% 
beyond that. It is not significant to the total central isovector component of the TPEP, since the 
contribution of family II is already fairly small [ffl]. The origin of such discrepancy still remains to 
be understood, and it might become relevant when more precision is asked for. 

For now on we ignore the above mentioned discrepancies and concentrate only on the effect 
of expanding (HB) or not (RB) the Lorentz-invariant loop integrals. In Fig.|| we show the ratio 
of the HB over the RB versions of the TPEP, projected into two specific partial waves. The 
illustrates the general behavior of partial waves with total isospin T = 1, which is similar to the 
LO+NLO terms in the HB expansion of the triangle integral (Fig.|l|). On the other hand, the partial 
waves with T = 0, represented by 3 Gs, are more sensitive to the HB expansion, with a discrepancy 
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of 10-30% between 2 and 5fm, and increasing to almost 50% at r = 15fm. This is due to the 
(isoscalar)— 3(isovector) structure in those waves, each term having nearly the same magnitude. 
Although in Fig.|3] we used the set EM 02 from Table [I[ these qualitative results do not change with 
the choice of LECs. 
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Figure 3: Ratios between HB and RB formulations of the TPER for G4 and G5 partial waves. 



LEC 


BM00 


Nij 03 


EM 02 


Set IV 


Set T] 


C\ 


-0.81 


-0.76 


-0.81 


-0.81 


-0.81 


C2 


8.43 


3.20 


3.28 


3.28 


3.28 


C3 


-4.70 


-4.78 


-3.40 


-3.20 


-3.80 


C4 


3.40 


3.96 


3.40 


5.40 


4.50 



LEC 


FMS 98 


di +d 2 


3.06 




-3.27 


d 5 


0.45 


di4 — d\5 


-5.65 



Table 1: Values for the nN LECs, from several authors: BM 00 [§], Nij 03 @, EM 02 0], Set IV [[TJ, and 
Set 77 [12\. The numerical values of the couplings di's, which have a small influence on NN observables, 
were taken from FMS 98 |J|]. 



4. Phase shifts and deuteron properties 

In Ref. [Q] phase shifts for some peripheral waves (F, G, and H) were computed using the 
available 71N LECs. The short-distance singularities were regulated with the same configuration- 
space cutoff used by the Argonne potentials [Q], [1 — exp(— cr 2 )] 4 xTPEP, with c = 2.0fm~ 2 . That 
work has shown that the differences between the RB and HB results for these waves are negligible. 
That happens due to the large superposition of the one-pion exchange potential (OPEP), making 
the difference in Fig|3] hard to observe. More pronounced is the dependence of the phase shifts 
with the set of LECs. Fig^ illustrates that with two selected partial waves. 



Results for lower waves and deuteron properties were calculated in Ref. [12] using the renor- 
malization method developed by the Granada group Jl6[]. The central idea of the method is that 
the behavior of the potential at short distances fixes uniquely the form of the wave functions at the 
origin 1 . A repulsive potential prevents the two particles to get close to each other, which causes 

1 The short-distance singularity of a potential is closely related to the divergent ultraviolet behavior of an interaction 
that is valid only at low energies. In order to have some predictive power, such effective theory has to be properly 
renormalized. 
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Figure 4: 1 F$ and 3 7*4 phase shifts, using LECs from Table 1. For comparison, we also plot results from the 
Nijmegen partial wave analysis [ |l5| | (open circles). 



the wave function to decrease exponentially in magnitude as r goes to zero. In the attractive case, 
on the other hand, physical conditions have to be imposed to avoid the particles to collapse at the 
origin with infinite velocity. In the Granada method, this is achieved by demanding the wave func- 
tion to satisfy specific boundary conditions at the origin. There the wave function is specified up 
to a phase, which is determined with an input of a physical quantity. For the generalized case of 
N coupled channels, a potential diverging at the origin as U(r) — ► m^Qn/r 11 , with C„ a matrix of 
generalized Van der Waals coefficients, is diagonalized via an unitary matrix G 



m N C„ = Gdiag(±/?!- 2 , 



,±R" N 2 )G 1 



(4.1) 



with Rj constants with length dimension. The plus (minus) sign corresponds to the case with a 
positive, attractive (negative, repulsive) eigenvalue. At short distances the solutions are given by 
u(r) — ► Gu, where u is a column vector with components u\ ,±(r), • • • ,UN,±(r). The attractive and 
repulsive cases behave respectively as 



Ui_(r) -> C,, — 



n/4 



sin 



- +<Pi 



H,+ \ 



n/4 



exp 



i-i' 



(4.2) 



(4.3) 



Here, <p,'s are arbitrary short distance phases which in general depend on the energy. There are as 
many short distance phases as short distance attractive eigenpotentials. Orthogonality of the wave 
functions at the origin yields the relation 



N 
i=\ 



r=0 



£C0S(<P;(£) -<?,(»), 
i=l 



(4.4) 



where A < is the number of the short distance attractive eigenpotentials. The <p,'s are fixed by 
low-energy (long-distance) input, which are the phase shifts close to threshold. The latter has a 
well-known low-energy expansion, the leading term being the inverse of the scattering "lengths". 
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For coupled channels with total (orbital) angular momentum j (j ± 1) the threshold behavior in the 
Stapp-Ypsilantis-Metropolis convention is 

8)U - & ' <, , 8)U - k^/d% , Sj -> k 1 ' ' aV . (4.5) 

Results for deuteron properties are given in Table |2[ For the RB potential using Set IV, the 
asymptotic D/S ratio r\ is a prediction, which overshoots the D-state contribution. We noted that 
this can be remedied by slightly changing the chiral couplings C3 and C4 to reproduce the experi- 
mental value of rj (Set 77). 



Set 


/(for 1 ) 






r<f(fm) 


Qd(fm 2 ) 


Pd 


OPE 


Input 


0.02634 


0.8681(1) 


1.9351(5) 


0.2762(1) 


7.88(1)% 


HB Set IV 


Input 


Input 


0.884(4) 


1.967(6) 


0.276(3) 


8(1)% 


RBE Set IV 


Input 


0.03198(3) 


0.8226(5) 


1.8526(10) 


0.3087(2) 


22.99(13) % 


RBE Set 17 


Input 


0.02566(1) 


0.88426(2) 


1.96776(1) 


0.2749(1) 


5.59(1) % 


Nijmll 


0.231605 


0.02521 


0.8845(8) 


1.9675 


0.2707 


5.635% 


Reid93 


0.231605 


0.02514 


0.8845(8) 


1.9686 


0.2703 


5.699% 


Exp. 


0.231605 


0.0256(4) 


0.8846(9) 


1.971(6) 


0.2859(3) 





Table 2: Deuteron properties, for several NN potentials. [|l2|. 



A surprising result of the RB potential with Set r\ is that the number of required low-energy 
inputs to renormalize the Schrodinger equation is significantly less than in the HB case, nearly a 
half (see Table II of Ref. Jl2|]). This is a consequence of the different behavior of the RB potential 
near the origin, V (r) ~ 1 /r 1 , which resembles a relativistic Wan der Waals force, in contrast with 
the typical 1 /r 6 of the HB counterpart. There is also see a clear improvement in the deuteron 
properties and on the energy dependence of several phase shifts, some of them shown in Fig. 
U There are, however, some partial waves ('So, x P>i, F3, and ^Pi) where agreement with the 
Nijmegen partial wave analysis gets slightly worse. Compared to the HB, the 'So, 3 Pi, £2, D3, and 
3 G5 RB results are not satisfactory, but in other cases they are either consistent or better, remarkably 
for [ D 2 , 3 P , 3 Pi, 3 £>i, and e { ph. 



5. Summary 

The program of constructing a nucleon-nucleon interaction based on the ideas of Becher and 
Leutwyler was presented. We pointed out the problem with the usual HB formalism in the one- 
nucleon sector and its consequences to the nucleon-nucleon system, namely, that it cannot account 
for the correct asymptotic description of the TPER If one performs a naive 1 /m# expansion of our 
results (which is not correct) one should, in principle, recover the HB expressions. In doing so, we 
noticed that there are three discrepant terms coming from two-loop dynamics. The origin of this 
discrepancy is not yet known and should be resolved if more precision calculations are demanded. 

Considering only the effects of the 1/m^ expansion in the covariant loop functions, we showed 
that the HB and RB versions of the TPEP can have considerable differences, in particular the partial 
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Figure 5: Selected partial wave phase shifts using Sets IV and 77, compared with the HB results (Set IV) 
and the Nijmegen partial wave analysis [OJ. 

waves with total isospin T = in the physically interesting region between 2 and 5fm. However, 
this effect is barely noticeable in peripheral phase shifts, due to the large one-pion exchange con- 
tribution. More pronounced is their sensitivity to the 7iN LECs, in particular, C3 and C4. 

To obtain the RB predictions for lower partial waves and deuteron properties, we used the 
renormalization method developed by the Granada group. We readjusted the LECs C3 and C4 to 
reproduce the asymptotic D/S ratio (Set T]), bringing other deuteron properties in agreement with 
experimental data and other potential model predictions. We also obtained an overall good agree- 
ment of the calculated phase shifts with the Nijmegen partial wave analysis, in a similar way as in 
the HB case. However, such agreement was obtained with significantly less counterterms (nearly a 
half) than the HB potential. That could be an advantage in determining the nucleon-nucleon LECs 
and bringing the RB potential into a form that could be eventually used in nuclear calculations. 
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